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We discuss a two-point particle irreducible (2PPI) approach to many-body physics which relies 
on a renormalization group (RG) flow equation for the associated effective action. In particular, 
the general structure and properties of this RG flow equation are analyzed in detail. Moreover, 
we discuss how our 2PPI RG approach relates to Density Functional Theory and argue that it can 
in principle be used to study ground-state properties of non-relativistic many-body systems from 
microscopic interactions, such as (heavy) nuclei. For illustration purposes, we use our formalism to 
compute the ground-state properties of two toy models. 



I. INTRODUCTION 

Originally, Density Functional Theory (DFT) has been 
invented by Hohenberg and Kohn to study atoms with 
many electrons in an efficient way [HIl]. Since then, DFT 
has indeed been successfully used to study a large variety 
of quantum systems with many degrees of freedom, rang- 
ing from electronic systems over ultracold Fermi gases to 
(heavy) nuclei. 

For heavy nuclei, DFT remains currently to be the 
only feasible approach for a calculation of ground-state 
properties [3]. In fact, the nuclear energy density func- 
tional approach represents a very active research field, 
also documented by the impressive efforts undertaken by 
the UNEDF SciDAC Collaboration [1|5]. 

The application of DFT to the nuclear many-body 
problem has been very successful in recent years. Apart 
from conceptional advances aiming at, e. g., ab-initio 
studies of heavy nuclei, nuclear DFT studies provided 
us with a universal understanding of properties of nuclei 
which is relevant for a large variety of applications, see, 
e. g., Refs. [BHS] and also Ref. [TU] for a review. 

For some time, DFT approaches have been based on 
fitting the parameters of a given ansatz for the density 
functional such that one reproduces the experimentally 
determined values of the ground-state properties of var- 
ious heavy nuclei |llj . The resulting density functionals 
have then been employed to describe ground-state prop- 
erties of other heavy nuclei. In recent years, there have 
been many attempts to give microscopic constraints on 
the nuclear energy density functional employing, e. g., 
(improved) density matrix expansions [TU [13] . The lat- 
ter approach has been further pursued and used to derive 
a (microscopic) nuclear energy density functional from 
chiral two- and three-nucleon interactions iT4l [T5] . On 
the other hand, the density-matrix expansion has been 
tested against ab initio calculations of trapped neutron 
drops 1Q\. These developments based on density-matrix 
expansions might be viewed as the beginning of a new era 
in the context of the nuclear energy density functional 



approach. Moreover, the future construction of density 
functionals will certainly benefit from various different 
approaches and complementary studies, ranging from a 
direct optimization of energy density functionals to stud- 
ies of the equation of state of nuclear matter, see, e. g., 
Refs. [TTHH]- At the present stage, however, our under- 
standing of the constraints for energy density function- 
als emerging from microscopic nuclear forces as well as 
the (direct) relation of the energy density functional to 
these forces is not yet fully complete and requires fur- 
ther research. An ab-initio renormalization group (RG) 
approach to DFT could complement and extend these 
efforts [211 [23] since it may open up the possibility to di- 
rectly compute ground-state properties of (heavy) nuclei 
from the underlying microscopic nucleon-nucleon interac- 
tions in a systematic fashion. In particular, such a func- 
tional RG approach seems to be promising since it allows 
to directly study the change of the energy density func- 
tional under 'RG transformations', e. g., from a weakly- 
interacting or even non-interacting system (starting point 
of the RG flow) to the fully interacting system, namely 
the nucleus under consideration [22 [23]. Moreover, such 
an RG approach is advantageous as it can be directly re- 
lated to the underlying path integral for which many sys- 
tematic approximation schemes are known, ranging from 
perturbative schemes to non-perturbative resummation 
techniques. For reviews and introductions to DFT ap- 
proaches using the path-integral formalism, we refer the 
reader to Refs. [Ml^ . 

From a very field-theoretical point of view, the ob- 
ject of interest is the so-called (quantum) effective action 
which can be derived from the (exponentiated) path inte- 
gral by introducing source terms and performing then a 
Legendre transformation with respect to theses sources, 
see also our discussion below. Such a construction is well- 
known in statistical physics where the Gibbs free energy 
of, e.g., a spin system, is computed from the underlying 
partition function by means of a Legendre transforma- 
tion. In quantum field theory, it is possible to work along 
these lines and couple the sources to the (physically) rel- 



evant "degrees of freedom" . For DFT, this means that 
we couple the sources to the densities associated with 
the quantum fields and then perform a Legendre trans- 
formation with respect to these sources. The effective ac- 
tion resulting form such a procedure is, strictly speaking, 
called a two-particle point-irreducible (2PPI) effective ac- 
tion. This 2PPI effective action is directly related to the 
energy density functional introduced by Hohenberg and 
Kohn, see our discussion below and also Refs. [m[3S] for 
a review. 

At this point, we have traced back the computation 
of the energy density functional to the computation of 
the path integral of the underlying theory defined by the 
microscopic interactions which itself represents an inher- 
ently difficult problem. Due to the relation of the en- 
ergy density functional and the path integral, however, 
we can utilize powerful existing tools for the computa- 
tion of the associated path integral. For example, we 
could employ ab-initio Monte-Carlo (MC) calculations 
or RG approaches. In the present work, we consider 
an RG approach to DFT that has been put forward in 
Refs. I2UI23J. For a more general discussion of the prop- 
erties of nPPI effective actions, we refer the reader to 
Ref. rSB]. 

In the present work, we give a detailed discussion of 
the RG approach to DFT introduced in Refs. [21 E5] . 
In Sect, [nj we present a general discussion of the struc- 
ture and properties of these DFT-RG fiows, including 
their connection to perturbation theory as well as the 
Hartree and Hartree-Fock approximation. In Sect. |IIl| we 
then apply our DFT-RG approach (2PPI-RG) approach 
to two simple toy models for which analytic solutions are 
known. In Sect. |IV[ we finally present our conclusions 
and outlook, including a concise discussion of our next 
steps towards an application our DFT-RG approach to 
the nuclear many-body problem. 



II. RENORMALIZATION GROUP APPROACH 
TO DENSITY FUNCTIONAL THEORY 

In this section, we give a more general discussion of 
field-theoretical aspects of the DFT-RG approach put 
forward in Ref. [221 and further discussed in Ref. (23. 



A. RG floAV equation 

In nuclear DFT, we are particularly interested in 
strongly-interacting many-body systems (far) away from 
the continuum limit where the ground-state density is 
inhomogeneousF] More generally speaking, we aim at 



-^ Note that strong interactions are by no means necessary to in- 
duce inhomogeneous ground states. For example, the ground- 
state density of A'^ non-interacting particles trapped in a har- 
monic oscillator potential is obviously inhomogeneous. 



a study of a finite system of fermions interacting via a 
non-local interaction which may be repulsive at short dis- 
tances and attractive at long range, as it is the case for 
the nuclear many-body problem and ultracold trapped 
Fermi gases. As discussed above, DFT has indeed proven 
to be useful for studies of these type of systems, see, e. g., 
Refs. [Tnil23[2H]- 

Before we discuss the DFT-RG flow equation and its 
properties, we give a brief summary of the underlying 
principles of DFT. To this end, we restrict ourselves to 
the following actionjj 

S[i^\^] ^ j j ^i(r, x) [a, - A + V{x)] ^,(r, x) 

+ \j j j i^l{T,x)i:l,{T\v')U{x,y)^„,{T\y')^„{T,x), 

where we have set 2m — 1 and introduced the short- 
hands /j. = /p dr^ Ix ~ I '^'^■^' ^ ^^ ^^^ space dimen- 
sion. The function V{x) denotes a (background) poten- 
tial. Moreover, we assume here and in the following that 
we sum over identical spin indices, if not stated otherwise. 
In the following we only consider two-body interactions. 
Higher-order iV-body interactions will be ignored but can 
be included straightforwardly in our DFT-RG approach. 
DFT is based on the famous Hohenberg- Kohn theo- 
rem [Tj. For a given interaction potential U, this theo- 
rem states that there exists a one-to-one correspondence 
between the ground-state density and the potential V{x) 
(up to an additive constant), at least for non-degenerate 
ground states. This implies that the ground-state den- 
sity rigs (uniquely) determines the ground-state wave- 
function of the A'^-body problem under consideration. 
The latter can therefore be considered as a functional 
of the density n. Moreover, the expectation value of any 
physical observable is determined by a unique functional 
of the ground-state density. In particular, this is true 
for the ground-state energy of the system and implies 
the existence of an energy density functional E[n\. The 
ground-state density rigs (x) can then be obtained by min- 
imizing E[n] with respect to the density: 

inf£;fnl. 



E, 



gs 



(1) 



Moreover, it can be shown that the energy density func- 
tional in the limit of vanishing external potential V^ the 
so-called Hohenberg-Kohn functional i?HKi is universal 
for a given interaction potential U: 



EuK[n] = E[n]- / d'^x 



n{x)V{x) . (2) 



These considerations can be generalized to the case of 
degenerate ground states. However, we leave aside a dis- 
cussion of the issue of F-represent ability, and also of N- 
representability, in the following. We refer the reader to, 
e. g., Refs. [29H3T] for a more detailed discussion. 



^ Throughout this work, we use the imaginary-time formalism. 
The extent of the imaginary-time axis can then be identified 
with the inverse temperature /3 = 1/T. 



The Hohenberg-Kohn theorem can be viewed as a start- 
ing point for an efficient description of many-body prob- 
lems. However, the theorem does by no means pro- 
vide a recipe for the computation of the Hohenberg-Kohn 
functionaL Similar to the IPI quantum effective action 
in conventional quantum field theory, the Hohenberg- 
Kohn functional consists of infinitely many terms. As 
we shall see below, the Hohenberg-Kohn functional is in- 
deed closely related to a specific type of effective action. 
This observation implies that it is in general not possi- 
ble to write down the exact Hohenberg-Kohn functional 
for a given many-body problem. Thus, an ansatz for the 
functional is required in order to determine the ground- 
state properties of a given many-body problem. Usually, 
it is difficult to find a systematic and stable approxi- 
mation scheme. For example, the simplest approxima- 
tion is the so-called local density approximation (LDA) 
which can be derived straightforwardly from the density 
dependence of the ground-state energy of the correspond- 
ing uniform many-body problem. The latter can, for 
instance, be computed with ab-initio MC simulations. 
It is indeed possible to show that LDA represents the 
lowest order in a systematic derivative expansion of the 
exact energy-density functional [32]. However, a low- 
order approximation of this type might only be justified 
in systems with weakly varying densities, such as ultra- 
cold Fermi gases with a large number of atoms in an 
isotropic trap |33j . For a general many-body problem, 
such a derivative expansion may have bad convergence 
properties. 

Now we would like to make contact between DFT and 
the effective action approach to quantum field theory, 
see , e. g., Refs. [MIES] for a more detailed introduction. 
The generating functional for our theory defined by the 
action S is given bjjj 



(3) 



^ qWI{j„}] 



In contrast to the conventional textbook approach to 
quantum field theory, namely the IPI formalism, we have 
coupled the external sources {Ja} to terms which are bi- 
linear in the fermion fieldsF] These bilinears play the role 
of composite bosonic degrees of freedom. 

Let us add a word on the issue of fixing the particle 
number in a path-integral approach. We can either intro- 
duce chemical potentials into the path integral to fix the 
numbers of the various particle species, say protons and 
neutrons, or we do not include chemical potentials but fix 
the particle numbers by choosing appropriate boundary 



^ Here, we have dropped an irrelevant normalization factor of the 
path integral. 

* Note that the formalism can be generalized by including sources 
coupled to, e. g., pairing densities. The latter might be consid- 
ered as convenient effective degrees of freedom to describe the 
ground-state properties of certain many-body systems. 



conditions for the equations of motion, as discussed in 
Ref. [34]. In a concrete DFT-RG study, we shall follow 
the latter approach to fix the particle number [33] . 

We now introduce the so-called classical fields p^ (t, x) 
which are defined as the (functional) derivative 
of W^[{>/ct}] with respect to the corresponding 
sources Ja-iT, x): 



Pa{T,x) 



5W[{J„}] 
5Ja (r, x) 



(4) 



Note that p^ is not only a function of t and x but also 
a functional of the sources {Ja\, i- e. p^ ~ Pa[{Ja}]- 
Clearly, these fields are related to the particle densities. 
In complete analogy to the textbook derivation of the IPI 
effective action, the 2PPI effective action is now defined 
as the Legendre transformation of W with respect to the 
sources J^'- 



nipA] 



sup 

{J.} 



-W[{J^}] 



Ja{T,x)paiT,x)} . (5) 



The so-defined 2PPI effective action r[{pcr}] determines 
completely the dynamics of the many-body problem and, 
up to a factor of /3, it can be associated with the energy- 
density functional mentioned above in the context of 
the conventional Hohenberg-Kohn DFT formalism. We 
add that the exact equivalent of the energy density 
functional as introduced by Hohenberg and Kohn can 
be derived similarly, if one introduces time-independent 
sources Jcr{x), see, e. g., Refs. [511 [35H55] . For a more 
general discussion on DFT in terms of a Legendre trans- 
formation, we refer the reader to Refs. |39» i40j . 

It is straightforward to show that the 2PPI effective 
action r[{po.}] does not dependent on the sources {Ja}'- 



ST[{Pa}] 
SJa 

On the other hand, we have 



0. 



5pa(T,x) 



= J a (t, x) 



(6) 



(7) 



Thus, the ground-state configuration {po-.gsl is deter- 
mined by this equation in the limit J<t — >■ OJj In other 
words, solving Eq. ([7| for the fields Pa{T, x) in this limit, 
we find the ground-state configuration {po-,gs}- We define 
the (time-independent) ground-state densities ngs,o-(x) as 
follows: 



^gs,a" y^} 



1 



dTPgs,o-{T,x) . 



(8) 



If the solutions {paiT, x)} turn out to be independent of 
the imaginary time r, then we have n^^gsix) = Pcr,ga{T, x). 



In the limit of vanishing sources, Eq. n\ represents the quantum 
equation of motion of the composite degrees of freedom pa ■ 



In our toy model studies to be discussed below, this is 
indeed the casel£] 

From the solutions {pa{T; x)}, we eventually obtain the 
ground-state energy Egs of the system: 



Egs ■■= lim -^mPa.gs}] ■ 



(9) 



This relation follows immediately from the spectral rep- 
resentation of the partition function Z: Z ^ ^^ q-I3e„ ^ 
and r ^ -W[J]\[j^^o} - \nZ[J]\y^^o}- From Eq. Q, 
we can also anticipate the relation between the effec- 
tive action F and the Hohenberg-Kohn functional -Ehk^ 

F ^ /3^HK. ^ 

At this point, we would like to add that the universality 
of the Hohenberg-Kohn functional i?HK follows from the 
fact that background potential can be absorbed into the 
source terms J^- by a simple shift, J^ ^^ Ja + Vy see 
Ref. [22]. Exploiting this observation, we find 

r[{Pa}] 

= ^KK[{Pa}]+Y.J dT jd''xV(x)p,{T,x), (10) 



where FHK[{/Ocr}] = F\/=o[{pct}]- We conclude that the 
functional FHK[{/3cr}] depends only on our choice for the 
interaction potential but not on the background poten- 
tial V . Recall that in our case the p^'s depend in gen- 
eral on the imaginary time r, in contrast to the stan- 
dard Hohenberg-Kohn functional which depends only on 
a time-independent density. In this respect, our present 
argument can be viewed as a trivial generalization of the 
original universality argument given by Hohenberg and 
Kohn. 

As discussed above, the computation of the effective 
action F (~ Hohenberg-Kohn functional) for a given the- 
ory can be an inherently difficult task. In the present 
work, we employ a DFT-RG approach that has been put 
forward in Refs. [551|2ni[3H]- For details on the derivation, 
we also refer the reader to these papers. For convenience, 
we only consider a system of N spinless fermions in the 
following. However, we stress that the derivation of the 
flow equation is by no means bound to such a theory 
but can be straightforwardly generalized to other non- 
relativistic theories. To be more specific, we consider a 
classical action of the following general form: 



Sx[r,i^] = 



r(T, X) [dr-A + Vx{x)] V'(t, x) + 



ij* (t, x)V'* (r', x')U{t, t'; x, xOV'It', x')tpiT, x) 



with 2m = 1. Here, we allow for a very general form of 
the two-body interaction. In most cases, however, one 
will restrict U to be of the form U ^ 6{t — t')U{x,x'). 
The parameter A £ [0, 1] denotes a dimensionless control 
parameter: For A = 0, the two-body interaction poten- 
tial U is turned off and we are left with an exactly soluble 
non-interacting problem, namely N fermions trapped in 
a given potential V\. For A = 1, the potential U is fully 
turned on and the potential Vx has assumed its physical 
form. For example, V\ can be chosen such that it plays 
the role of the trap potential in an experiment with ultra- 
cold Fermi gases. In studies of ground-state properties 
of self-bound many-body problems, such as nuclei, one 
chooses V\=i{x) = 0. Apart from the physical constraint 
at A = 1, the form of Vx is at our disposal and can be cho- 
sen such that, e. g., the initial non-interacting problem 
is simple to solve. For example, we could choose Vx=o to 
be a harmonic potential. 

We add that, from a field-theoretical point of view, 
the one-body potential Vx acts as a regulator function. In 
fact, the length scale associated with the potential Vx sets 
a momentum scale in the theory which removes infrared 
divergences in the appearing loop diagrams. In this sense. 



® In fact, this T-independence of pgs follows from our flow equa- 
tion | |11[ | to be discussed below, provided that we only consider 
an interaction potential U with a imaginary-time dependence of 
the form U ~ S(t — r'). 



the introduction of the parameter A renders the theory 
scale-dependent. In particular for the case Vx=i{x) = 
0, we shall also assume that the Fourier transform of 
the interaction potential U falls off sufficiently rapidly 
for large momenta to avoid the occurrence of ultraviolet 
divergences. 

Following Refs. [22l[23], it is now straightforward to de- 
rive the RG flow equation for the 'scale-dependent' 2PPI 
effective action Fa [p] ■ By taking the derivative of F^ [p] 
with respect to A, we flnd 



dxrx[p]^{dxVx)-p+lp-U-p 



-TtU ■ 
2 



SpSp 



(11) 



(12) 



where the dot represents a shorthand for 

AB= I dr I d'^xA{T,x)B(T,x) 
and the trace Tr stands for 



TrM(T',f' ,r,f) = / dr / d'^xM{T,x,T,x) . (13) 



This functional differential equation (11) describes the 



flow from the non-interacting system defined at A = to 
the interacting theory defined at the physical point A = 1. 
The effective action Fa=o[p] associated with an exactly 



soluble (non-interacting) TV-body problem determines 
the initial condition of the RG flow. 

In the terminology of many-body physics, the sec- 
ond term in Eq. (II) can be identified as the so-called 



Hartree term. The third term on the right-hand side de- 
pends on the scale-dependent density-density correlator 
8'^T x[p] / [5 p5 p) and includes all other corrections to the 
effective action, also so-called Fock contributions, see also 
our discussion in Sect. Ill Dl 



We observe that the flow equation (II) has a sim- 



ple one-loop structure as it is the case for the RG flow 
equation for the IPI effective action derived by Wet- 
terich [4T1|42]. However, this does not mean that we can 
only capture one-loop corrections with this flow equation. 
On the contrary, by solving the functional differential 



equation ( II ), we automatically include corrections of ar- 



bitrarily high-order, see also our discussion in Sect. II C 



In fact, no approximations are involved in the derivation 
of the flow equation (II). In particular, the derivation 



of this flow equation does not require that the interac- 
tion strength is small. The equation is exact, provided 
we only allow for two-body interactions in the classical 
action S'aQ 

For a general many-body problem, it is not possible 
to solve the flow equation (II) for the 2PPI effective ac- 



tion (corresponding to the Hohenherg-Kohn functional) 
exactly. Thus, we need to find (systematic) approxi- 
mation/truncation schemes. Examples for such approxi- 
mation schemes are the gradient expansion of the effec- 
tive action or an expansion of Tx[p] about the ground 
state Pgs,A: 



Y\[P] =rA[pgs,A] + 7, 



{p{t, x) - PgsAr, x)) rf ^ [Pgs,\]iT, x; r', x') {p{t' , x') - Pgs,A(T', x')) -h . . . . (14) 



Plugging this expansion into the general flow equa- 
tion ( II ), we obtain the flow equations for rA[pgs,A]j Pgs.x, 



and T^^'[pgs.\] with n > 2. We then find that the flow 
of Pgs.A depends on Pgs,x itself but also on r\ [pgs.x] 



and rf^[pgs,A]- 
point function T^ [p^ 



ri%: 



gs,Ajj 



.(4) 



The flow equation for the two- 
'gs,A] depends on p^^^x, ^x iPgB,x]: 



and P^ [/Ogs,A]- In general, we find that 



^(n)r 



-i(m) r 



the flow of P^' [Pgs.A] depends on Pgs,x and P}^ [/3gs,A] 
with m = 2, . . . ,n,n + l,n + 2. 

In the terminology of quantum field theory, the expan- 
sion ( 14 ) corresponds to a vertex expansion and will be 



discussed in detail in our toy model studies in Sect. |III| 
Loosely speaking, the n-point correlation function is re- 
lated to the expectation value of n density operators: 



^(")r 



,a](ti,^1; • ■ • \Tn,Xn) 

~ {p{'ri,xi)p{T2,X2)- 



■p(T„,f„))gs, (15) 



where 



p(t,x) =-0*(T,f)V'(r,f). 



(16) 



For brevity, we have dropped the subtraction of the 
disconnected contributions on the right-hand side of 
Eq. (15), see Eq. (19 1 for a more rigorous expression 
for n = 2. At this point, we would like to stress that the 



vertex expansion given in Eq. ( 14 ) is an exact expansion 



of the functional Pa [p] about the ground state and should 



^ A generalization of this flow equation to include the effects of 
higher n-body operators is straightforward. In principle, this 
requires to also multiply the 3—, 4—, . . . body-interaction terms 
with the control parameter A. 



I 

by no means be confused with the local density approx- 
imation. Recall that the expansion coefficients, namely 
the n-point functions, are functions of the time-like and 
spatial coordinates. 

Before we discuss the features of our RG flow equation, 
we would like to emphasize that our present approach al- 
lows us to directly compute the energy density functional 
from the microscopic interactions of the theory. To be 
concise, it opens up the possibility to derive the energy 
density functional (and the ground-state properties) of, 
say, heavy nuclei from chiral effective field theory inter- 
actions [43*-H7] as the latter can be used to determine the 
interaction potential U in our flow equation. 



B. Excited States 

Up to this point, we have only discussed how to com- 
pute the ground-state energy. However, it is also possible 
to extract the energy of the excited states of the theory 
from the 2PPI effective action T[{p„}] = Pa=i[{p^}]. To 
this end, it comes to our rescue that we consider the sys- 
tem at a finite temperature T. Using now the spectral 
representation of the partition function Z, we find 



r[{p.,gs}] = -in E' 



'I3E„ 



const. 



(17) 



In order to obtain the energy of the excited states En 
(n > 0, Egs = Eq), we could simply compute P[{P(T,gs}] 
as a function of f3 and then fit the result to the func- 
tional form given on the right-hand side of Eq. (17 1 with 



the energies of the excited states En as the free fit pa- 
rameters. For sufficiently low temperatures (i. e. suffi- 
ciently large /3), it is reasonable to assume that only the 



lowest- lying states are occupied. In a fit, we would then 
only include a small finite number of fit parameters cor- 
responding to these states. However, such a procedure 
may still turn out to be impractical from a numerical 
point of view. 

In order to extract the energy of the first excited state, 
we do not require a fit procedure involving the spectral 
representation of the partition function. In fact, the first 
excited state Ei can be extracted from the two-density 
correlation function evaluated at the ground state. 



.(2) 



rr;,(r,f;r',f') 






{P<J,Es} 



(2) 

First, we note that F^ ^, is directly related to the (field- 
dependent) propagator Ga-,c7'- 



G. 



,.,(T,x;r',f') = (r5'(^,2f;r',f')) ' 



(18) 



On the other hand, the propagator can be written in 
terms of expectation values of density operators: 

Ga-,a-'\'''i X;T ,X ) 

^{p„{t,x)p„>{t' ,x'))~'{p„{t,x)) {p„,(t',x')) . (19) 

The quantity pa- is the straightforward generalization of p 
defined in Eq. (16 1. From this expression, it follows that 



(El - E, 



) — — lim 

/3-i.oo 



lnG,,,(0,x;/3,f'), (20) 



independent of our choice for the space-time coordinate- 
pairs (t, a?) and (t',x'). Thus, we can extract the en- 
ergy of the first excited state from the knowledge of the 
two-point correlation function once we have computed 
the ground-state energy using Eq. ([9]). In principle, it 
is also possible to project on higher-lying excited states 
by considering the expectation values of properly chosen 
operators. In any case, our discussion shows that excited 
states are also accessible within our 2PPI effective action 
approach. 



C. Perturbation Theory 

Let us now discuss the properties of our DFT-RG ap- 
proach in more detail. For convenience, we now restrict 
ourselves again to the case of spinless fermions, i. e. we 
shall consider the flow equation (11) in the following. 



However, our line of arguments also holds for any type of 
fields: fermions with spin or even scalar fields as we shall 
see in our toy model studies in Sect. |III| 

In order to recover perturbation theory from the flow 



equation ( 11 ), we first introduce a counting parameter uq 
into the t 



where 



leory as follows: 

C/(t, t'; X, x') = UoU{t, t'; x, x') , (21) 



U{t, t'; X, x') = MS, x')5{t - r') . (22) 



Here, the full dependence on the space coordinates has 
been absorbed into the function 4)\j. For convenience, 
we assume that the dependence on the imaginary time is 
given by (5(r — r') and that 0u is a dimensionless func- 
tion. Now we can introduce a dimensionless counting 
parameter uqj 



Mo = ii/Uo . 



(23) 



Here, iy denotes a typical length scale associated with 
the background potential V. For example, we have 
Hv = 1/%/^: if we choose a harmonic background po- 
tential: V{x) — (l/2)w^x^. Using the length scale ^y, 
we can also introduce the dimensionless ground-state en- 
ergy, density, as well as dimensionless correlation func- 



tions F 



and 



("). 



E, 



gs 



^vEg 



Pgs — (-vPg. 



p(n) _ ^2np(n) 



(24) 



(25) 



where F^"^ — (5"F/(5p". Note that (.y'^P is dimensionless. 
The dimensionless ground-state energy, density, and 
correlation functions can be expanded in powers of uq: 



Egs — 2_^ ^nUQ , 
n=0 

oo 

Pgs = / ^ ^n^o ' 
11=0 

oo 



(26) 

(27) 
(28) 



n=0 



..M 



The coefficients eo, t'o, and 7g are determined by the 
initial condition of the flow equation, i. e. by the func- 
tional F>,=o of the non-interacting theory. At A = 0, all 
the expansion coefficients with n > are identical to 

zero. However, all the coefficients e„, i^n, and 7„ with 
n > depend implicitly on A and may therefore be gener- 
ated dynamically by quantum corrections. Note that the 
coefficients Vn also depend on t and x. Analogously, the 
coefficients 7^ depend on m pairs (r, x) of space-time 
coordinates. 

Let us now analyze the perturbative expansion of the 
ground-state energy. To this end, we apply the expan- 
sion ( [l4| ) about the current ground-state to our general 
flow equation (11). This yields the following equation 
for FA[pg 



i^gsj 



^aF 



A-L A[ 



.AJ = 



(dxVx) ■ Pgs,\ + ^Pgs,\ ■ U ■ pgs,\ 



-TtU ■ 



s'r^ip] 



SpSp 



(29) 



Pgs 



Plugging the expansions ( 26 )-( 28 ) into this flow equation 
and noting that -Egs, a = ^[Pgs,x\/l3 for (3 —>■ 00, we can 
derive flow equations for the coefficients e„ by simply 



comparing the left-hand side and right-hand side order 
by order in our expansion in powers of uq . 

Let us now distinguish between two cases, namely the 
case with a A-dependent background potential and the 
one with a A-independent background potential. We be- 
gin our analysis with the latter case. The first term on 
the right-hand side of Eq. ( 29 ) then vanishes identically. 



To the flow of the coefficient eoj no term on the right-hand 
side of Eq. ( 29 ) contributes since they depend explicitly 
on the interaction potential. Thus, these terms can only 
contribute to the flow of the coefficients e„ with n > 0. 

In order to obtain the correct result for the ground- 
state energy £^gs (^ r[pgs]) in leading order in uq, we 



„(2) 



deduce from Eq. ( 11 ) that we need to compute the coef- 



ficient 7g which is associated with the term independent 



of Uq in the expansion ( 28 ) : 



TtU ■ 



S^Tx[p] 



Sp5p 



uqTiU ■ 7 



.(2) 



■Oiul) 



Interaction-induced corrections of the propagator do not 
contribute to the leading-order correction of the ground- 
state energy. This follows immediately from the fact that 
the third term on the right-hand side of Eq. ( 29 ) depends 
explicitly on the interaction potential U and therefore 
on Uq. Since the second term on the right-hand side has 
an explicit [/-dependence as well, the coefficient vq is also 
required to recover the result from perturbation theory 
at leading order: 



Pgs.X ■ U ■ Pgs,A ~ Uq[vq-U ■ 1^0 j + 0{Uq) . 

These considerations can be continued successively. In 
general, we find that we need to compute pgs and the two- 
point function up to order Uq~ , in order to obtain the 
correct result for the ground-state energy up to order Uq . 

Let us now turn to the case with a A-dependent back- 
ground potential. We can follow the same line of argu- 
ments as in the previous case. Due to the presence of the 
term (dxVx) ■ Pgs,x, however, we now find that we need to 
compute Pgs,A and the two-point function up to order Uq 
in order to obtain the correct result for the ground-state 
energy up to order UqIj 

In summary, we have seen that the computation of 
the ground-state energy up to a given order uq requires 
that we also compute pgs and the two-point function up 
to a certain order. The required expansion order for the 



* In fact, this case is more subtle. Whereas it is apparent that 
we have to compute Pgs up to order Ug in order to obtain the 
correct result for the ground-state energy to the same order, one 
might naively expect that we only need to compute the two-point 
function up to order Mq~ , as in the previous case. However, the 
computation of pgs up to order Mq requires that we know the 
two-point function up to order Uq , in the case of a A-dependent 
background potential. Overall, we therefore need to compute pgs 
and the two-point function up to order Uq in order to obtain the 
correct result for the ground-state energy up to this order. 



latter depends on whether the background potential is A- 
dependent or not. 

As already mentioned above, the flow equations for 
the density and the two-point function (and therewith 
for their expansion coefflcients) can be obtained by ex- 
panding the general flow equation ( 11 ) about the current 
ground state Pgs,A and then projecting it on the corre- 
sponding quantities, namely rA[pgs,A]i Pgs,\i and the n- 



(")r 



point correlation functions T\ [pgs.x], see also Eq. (14 1 



We shall discuss this procedure in great detail in our toy 
model studies in Sect. |III| At this point, however, we al- 
ready would like to emphasize that the vertex expansion 
should by no means be confused with the perturbative 
series expansion discussed above. The associated expan- 
sion coefhcients are inherently non-perturbative quanti- 
ties. Below, we use the vertex expansion about the cur- 
rent ground-state since it is systematic and allows us to 
extract the RG equations for rA[pgs,A]j Pgs,Xi and the n- 
point correlation functions F^^" [pgs.A] in a simple manner. 



D. Hatree approximation 

The so-called Hartree approximation can be obtained 
from the general flow equation ( |11[ ) by dropping the third 
term on the right-hand side. The latter includes, e. g., 
the so-called Fock term. Dropping the third term on the 
right-hand side of Eq. ( 11 ), we can solve the flow equation 
for Fa[p] analytically. We find 



TaM = (Vx 



Va=o) ■ P+ ^P-U ■ p 



rx=o[p] , (30) 



where Fa=o [p] is simply the 2PPI effective action of the 
non-interacting (initial) system at A = 0. From the defi- 
nition of the ground-state, 



STxip] 



6p 



= 0, 



we obtain the implicit equation 



Pgs.A 



>-' 



AV^, 



srx=o[p] 



6p 



Pes, A. 



(31) 



(32) 



with AVx = {Vx — Vx=q)- Expanding the initial effec- 
tive action Fa=o[p] about the initial ground-state pgs,o = 
Pgs,A=o only up to second order, we find the following 
solution for Pgs.x'- 



Pgs,A - 



A, 



where F, 



(2) 
gs,0 



-u 



.(2) 



.(2) 



.(2) 
-gs,( 



AVx-F-Q ■ pgsfi 



(33) 



^gs,A=o[/^gs,A=o]- However, such a low- 
order expasion in p can only be meaningful if the in- 
teracting ground-state Pgs,A=i is close to the initial non- 
interacting ground-state Pgs,A=o- In general, it is difficult 



to judge a priori whether this is the casejj From a prac- 
tical point of view, we therefore have to include higher 
orders in the expansion of r;^^o [p] ^^^ analyze the con- 
vergence of the physical observables as a function of the 
expansion order. 

From our discussion, it follows that the Hartree ap- 
proximation already yields arbitrarily high orders in an 
expansion in powers of uq, both for the ground-state en- 
ergy and density. Recall that U ^ uq and E^s ^ ^[Pgs]- 
However, we would like to point out a shortcoming of 
the Hartree approximation which becomes apparent from 
our analysis. Taking into account our findings from 
Sect. |IIC[ we conclude that the Hartree approximation 
necessarily fails to reproduce the perturbative result for 
the ground-state energy, even at leading order. This is 
simply due to the fact that the third term on the right- 
hand side of Eq. (11) is missing in this approximation. 



This term depends explicitly on the interaction potential. 
Since S^T\/{SpSp) is in general not identical to zero, even 
in the non- interacting limit, this term generates terms 
which already contribute to the leading order in a per- 
turbative expansion. 



III. CASE STUDIES 

In this work, we refrain from an explicit study of 
ground-state properties of nuclei but rather present toy 
model studies to explain the theoretical formalism de- 
tailed in the previous section. The (classical) actions S 
underlying our toy model studies resemble the action un- 
derlying non-relativistic many-body problems in a few 
aspects, at least from a purely field-theoretical point of 
view. However, we also would like to emphasize that such 
toy model studies only represent a first step. Studies 
of simple self-bound many-body models, which are also 
much closer to nuclear physics from a phenomenological 
point of view, will be presented elsewhere [35] . 



A. Zero-dimensional Toy Model 

The simplest example for the application of our RG ap- 
proach is the computation of 'ordinary' integrals which, 
loosely speaking, corresponds to zero-dimensional field 
theory. In this case, the partition function correspond- 
ing to Eq. pi is an 'ordinary' integral of the forrrrj 



Z[J] 



dipe 



^S[^]+Ji,^ 



(34) 



^ Strictly speaking, the notion 'close' requires the definition of a 
measure on the space defined by the functions p. We shall skip 
this issue here. 
^^ Here and in the following, we drop again irrelevant normaliza- 
tion factors of the partition function. Moreover, we note that the 
quantities Z[J], r[p], . . . are no functionals but 'ordinary' func- 
tions in the present case. Nevertheless, we stick to our notation 
introduced in the previous section. 
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Figure 1. Effective action T[p] of the zero-dimensional toy 
model for the non-interacting case uo = and for ito = 1, as 
obtained from a direct calculation of the partition function. 
We have normalized the effective action such that r[pgs] = 
for Uo = 0. 



where 



sm 



2„;,2 



-W^V 



Mo 

24 



r 



(35) 



Here, J and ip are real-valued numbers rather than fields. 
Thus, the derivative terms in the action vanish identi- 
cally. Moreover, we have chosen U — uq/12 for the in- 
teraction potential and V\ = (l/2)(u;^ for the background 
potential. We add that zero-dimensional models have 
been already successfully employed to benchmark other 
field-theoretical methods, see, e. g., Ref. [48] . 

For J = J/o;^ < —1/2, the partition function can be 
given in closed form: 



Z[J]' 



V3-6JK1 



/ 3(1-2./)^ \ 
V 4fio J 



3(l-2J)^ 

e ^"'0 



,WIJ] 



/Uo 



(36) 



where uq = uq/uj'^ and K^ is the modified Bessel func- 
tion of the second kind of order z/, see, e. g., Ref. 1491 . 
For J > —1/2, the integral Z[J] can still be computed 
numerically. 

Using Eq. Q and taking the limit J — >■ 0, we can now 
compute the ground-state pgg = (';/'^)gs: 



Ps 



Pgs 



UJ UqKi 



2 3 



1- 



{uo + 3)Ki 



3 

4uo 



3 

Auo 
3 

This expression can be expanded in powers of uq'- 

Uo 
2 



2 " 



2^2 



11^3 



3 

4uo 



Oi 



^0, 



(37) 



(38) 



The ground-state 'energy' can be obtained directly from 
E^, = -{\nZ[0]-\nZ[0]l^^^) , (39) 
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Figure 2. Ground-state 'energy' E^s and ground-state 'density' pgs of the zero-dimensional toy model as a function of the 
(dimensionless) coupling uq. For comparison, we also show the results from a small coupling expansion at leading order (LO) 
and next-to-leading order (NLO). 



where we have normalized Egg such that it is zero 
for uo — > 0. For small uq, we find 






12 



11S| 
96 



0{ut) 



(40) 



Alternatively, we may compute the effective action T[p\ 
using Eq. ([5]) . From the minimization of the effective ac- 
tion, we then obtain pgg and Egg. Note that the computa- 
tion of the effective action requires that we solve p = p[J], 
as defined by Eq. ^, for the source J. The solu- 
tion J = J[p] needs to be plugged into the definition ^ 
of the effective action. For the case uq — 0, for instance, 
the effective action can be computed analytically. We 
find 

Tu,=o[p] - I (p^' - In i^TTpw^) - 1) . (41) 

Apparently, any global analytic ansatz for F^jq^o is bound 
to fail. However, a Taylor expansion about the ground- 
state PgsW^ = 1 (for Uq = 0) is possible and meaningful. 
For Uq > 0, the computation of the effective action can 
only be performed numerically. We find that the non- 
analyticity at p = persists in this case. In Fig. [T] we 
show the effective action T[p] as a function of p for uq — 1. 
We have normalized F[p] such that r[/9gs] = for uq — 0. 
For Uq = 1, the ground state pgs is found at pgs ~ 
0.750 . . . , in agreement with our analytic result (37). For 
the ground-state 'energy', we find ii^gs = 0.084 . . . , which 
also agrees with the analytically found value as obtained 
from Eqs. (36 1 and (39). We can compare these results 



with those from the effective action in the Hartree ap- 
proximation as derived from Eq. (30): 



"0,_ 2n2 ^ r , 1 



rHartrccW - ^ (P^ ) + ^ ^.=M + ^ ln(2^) . (42) 

The normalization has been again chosen such 
that FHartrcc[Pgs,Hartrcc] = for Uq -^ 0. For the grouud- 



state, we find 



/3gs,Hartroc(uo) = "q ^ (\/6uo + 9 - 3) 



1 



^-ii-i3+^«)-<«> 



The ground-state energy is then given by 



^gs("o) 



Hartree [Pgs, Hai'ti'coJ 



24 



"0 

144 



2592 



0\ 



(44) 



Clearly, the Hartree approximation does not reproduce 
the exact results, even in the small-coupling limit. In 
fact, it significantly underestimates the exact results 
for ii^gs and overestimates those for the ground-state pgg. 
In Fig. [21 we show our results for Eg^ and pg^. The com- 



parison ofthe exact results with the results from a small- 
coupling expansion at leading and next-to-leading order 
shows that the latter are only meaningful for -uq < 0.2. 
The results from the Hartree approximation are not in 
agreement with the exact results, neither at small cou- 
pling nor at strong coupling. For uq = \. for example, 
the relative error of the Hartree expansion amounts to 
about 57%. 

For illustration purposes, we now compute Eg^ and pgs 
with our flow equation (11). In the present case, it as- 
sumes a simple form: 



d\^\{p\ = ^^^0^^ 



5p6p 



(45) 



In order to solve this equation, we expand V[p\ about the 
current ground-state Pgs,\^ see also Eq. (14): 



Y\[P] =rA[pgs,A] + 



E 

n=2 






[Pgs,A](p-Pgs,A)",(46) 
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Figure 3. Effective action r[p] = rA=i[p] of the zero- 
dimensional toy model as obtained from our RG approach 
for various values of iVmax for fixed mo = 1- Again, we 
have normalized the effective action such that rA=i[pgs] = 
for uo = 0. 



Figure 4. _Egs as a function of uo as obtained from our RG 
study. Note that our results for A^max = 4 are already almost 
indistinguishable from the exact values for i?ga on the scale of 
the plot. 



Note 

that ^Ar^^ [pgs.A] = and r^^^gbgs.A] = by definition. 
Plugging this expansion into the flow equation ( 45 ) then 
yields a tower of flow equations for i?gs,A ~ ^\{P^&,>\^ 

The initial 



^(n) r 



Pgs.A: and the n-point functions r^"'^[pgs a] 

conditions for these flow equations can be extracted by 

expanding FflQ^ob] about its ground state Pgs. 

For A'niax = 2, for example, we find the foUowing set 
of coupled ordinary first-order differential equations: 



d\Egs,\ — 



1 - 4 

24"°^ 


pi + (rf[p.s,A])"' 


,(47) 


- — U0'^Vgs(ri^^bgs,A]) , 


(48) 


1 - 4 

12"°" • 




(49) 



^APgs.A = 



Note that -/Vmax = 2 is sufficient to exactly repro- 
duce the perturbative results for the ground-state en- 
ergy Sgs at leading order. In order to correctly repro- 
duce the leading order of the perturbative expansion 
of the ground state pgs and the n-point functions of 
higher order, we need to increase the truncation order 
beyond iVmax = 2. From our general discussion above, it 
follows that TVinax = 4 is sufficient to correctly reproduce 

(2) 

the perturbative series of pgs and r\ [pgs, a] at leading or- 
der. Moreover, we also recover the correct results for Egg 
at next-to-leading order with iVmax = 4. 

We would like to add that the set of fiow equations 
for TVmax = 2 can be solved analytically. We refrain from 
giving the explicit result here. We only state that the 
ground-state energy behaves as 



Egsiua) -- InuQ. 



(50) 



This is in accordance with the asymptotic behavior of the 
exact result, only that the coefficient of this term is not 



reproduced correctly for A'max = 2. 

In Fig. |3J we present our results for T[p] = Fa=i[p] 
as a function of p for various values of A'max for 
fixed uq = 1. Note that the radius of convergence of 
our expansion (46 1 about the ground state is finite. 
In fact, we find Tp/uj'^ = 1 for the (dimensionless) ra- 
dius of convergence in the case uq = 0. Our numeri- 
cal results for finite uq are in accordance with this re- 
sult. In fact, we do not observe a convergent behavior 
around pcu'^ « 2 for increasing A^max- On the other hand, 
our results for F[p] nicely approach the exact results for 
\p — Pgs\^'^ ^ 1 when A'max is increased. In order to com- 
pute T[p] for yow^ > 2, we could employ Taylor expansions 
around various different points with overlapping regions 
of convergence. This would be of importance, for exam- 
ple, when we expect that T[p] develops various minima. 

In Fig. Ill we present our RG results for Egs as a func- 
tion of Mq for various values for A^max- We find that our 
results for A'^max = 4 are already in very good agree- 



ment with the exact results for 



uo 



< 



1 and that we ap- 



proach the exact results from above for increasing A'max- 
For large values of the coupling uq , we observe that we 
need to go to higher truncation orders in order to repro- 
duce the exact results. Assuming that we do not know 
the exact solution for a given value of uq, these find- 
ings imply that we have to compute Egg as a function 
of A^max and check numerically the convergence of this 
function. This is illustrated for uq = f in Fig. [5] The 
solid (blue) line shows the result from a fit of the RG data 
for A^inax ~ 2, 3, 4, 5, 6 to the (empirical) three-parameter 
ansatz 



^gT(^n 



E. 



fit 



aiC 



-Ol2N„ 



(51) 



Here, E^^, ai, and a2 are the three fit parameters. We 
obtain £'g* w 0.0833 which is about 2% smaller than the 
exact result E^s ~ 0.0846. 
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Figure 5. i?(Sgs) = E^^/E^^'"'^ as a function of iV^ax 
for Wo = 1. The solid (blue) line represents a fit of the RG 
data to the (empirical) ansatz (51 1, see main text for details. 



B. 



Quantum Anharmonic Oscillator: 
One-dimensional Toy Model 



Up to this point, we have only discussed a zero- 
dimensional model which essentially boils down to or- 
dinary calculus. Let us now discuss a simple one- 
dimensional field-theoretical toy model described by the 
following classical action: 



S 



1 



P 



2 JO 

1 rf^ 



dTi;{T)[-dl+CJ^]2P{T) 



+ - dr dT'^{T)xP{T')U[T~T')i![T')ij{T) , (52) 





where iP{t) is a real- valued field and 

U{t-t') = ^Uo5{t-t' 

From a quantum mechanical point of view, this classical 
action describes nothing but the quantum anharmonic 
oscillatorpM Note that V'(''") can be associated with a 
time-dependent coordinate. We have chosen the normal- 
ization factor of the interaction potential U such that it 
corresponds to the standard choice for this model, see, 
e. g., Ref. HO]. 

Much is known about this model. In fact, it is straight- 
forward to diagonalize the Hamilton operator of this 
model numerically without any approximation, see, e. g., 
Ref. 51j. Thus, we have again an exact solution at hand 
which allows us to benchmark our RG results. We would 
like to add that this model has already been employed to 
benchmark other RG approaches, such as the IPI RG 
approach [S2J [53] and also 2PI approaches ^A\. The 



large-coupling limit of this model has been studied us- 
ing so-called large- iV techniques, see, e. g., Ref. [50j- In 
the latter approach, it is straightforward to show that 
the ground-state energy scales as 

-Bgs - uul (54) 

for wq 3> 1. Here, uq is the dimensionless coupling con- 
stant which is defined af^ 

Wo 
,,,3 • 



Uo 



(55) 



In the following we use our 2PPI RG approach (DFT- 
RG approach) to study the ground-state properties of 
this model. To this end, it is convenient to expand the 

ground-state Pgs,x{T) and the n-point function P^" [Pgs,A] 
in a Fourier series. For the ground-state, we choose 

1 " 



Pgs,A(r) = 



/3 



E pi"^e" 



(56) 



where a;„ = 27rn//3. Note that the fields "ipir) are real- 
valued fields and therefore obey periodic boundary con- 
ditions in the imaginary-time direction. 

(2) 

Rather than working with the two-point function P^ , 
we shall use the propagator Gx{t,t') from now on. The 
latter is defined as the inverse of the two-point function: 



/3 



.(2 



dT"GxiT,T") P^^^[pg,,;,](r",T') =<5(r-r').(57) 



Note that G\ should not be confused with a single- 
particle propagator. However, both propagators are re- 
lated in simple terms, see our detailed discussion in Ap- 
pendix IX] The propagator can again be expanded in a 
Fourier series: 



1 



Gx(t,t') = - 



E ^E 



qM^U 



{m,n) 



(53~> The matrix G\ ' is diagonal, i. e 



G 



{niju) 



(3G[" 



(58) 



(59) 



Higher n-point functions can be expanded accordingly: 

Pi"^[Pgs,A](ri,...,r„) 



TTIti— - 



A 



iEr=i' 



(60) 



In order to derive the flow equations for the Fourier 
coefficients, we again employ our vertex expansion, see 
Eq. ( 14 ) . Plugging the latter into the general flow equa- 
tion ( 11 ), we eventually obtain the flow equations for the 



Fourier coefficients. The initial conditions for these RG 
equations are determined by the non-interacting problem 
{uq — 0), see Appendix [A] 



^^ From a field-theoretical point of view, it corresponds to a so- 
called i/>^-theory in one dimension. 



^^ Note that the dimension of the coupling in units of uj difli'ers from 
the one discussed in Eq. | |23| l. This can essentially be traced back 
to the fact that the time derivative appears quadraticaUy in the 
action (|52| of our present model rather than linearly. 
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Figure 6. Comparison of the ground-state energy Eg^ as a 
function of uo as obtained from various different approaches. 
The RG results are compared with the exact results and the 
perturbative results at leading order (LO) and next-to-leading 
order (NLO). The analytic solution (681 for the leading-order 



RG approximation (RG, F^'^') overestimates significantly the 
ground-state energy. For details on the RG approximations 
of higher order (RGs, F^'^'; RGi, F'*'), we refer the reader to 
the main text . 



1. Leading- order approximation 

Let us begin with a discussion of the lowest non-trivial 
approximation which is given by dropping T^ [pgs,\] as 
well as their RG flows for n > 3, i. e. we consider the 
case A'lnax — 2. For the ground-state energy Egs,\, we 
then find 



d\Egs,\ — 



24 



1 -(0) 



oo 

E 

m — — oo 



1 



G 



(m) 



(61) 



Here, we have introduced the following dimensionless 
quantities: 



-(m) 2 (™) 



(62) 
(63) 
(64) 
(65) 



In terms of these quantities, the flow equations for the 
Fourier coefficients associated with the ground state pgg 
read 



CaPa 
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Finally, the RG flow of the propagator is determined by 
the following set of equations: 



SaGI"^ 



1 

12 



uoiG 



^(m) 



(67) 



Clearly, Eqs. (61 1, (66) and (67) represent an infinite set 



of fiow equations. For the present leading-order approx- 



Figure 7. Comparison of the ground-state pgs as a func- 
tion of Uo as obtained from various different approaches. The 
analytic solution ( 72 1 for the leading-order RG approxima- 



tion (RG, r"') is not in accordance with the exact results, 
neither in the small coupling-limit nor for large values of the 
coupling, see main text for details on the higher-order RG 
approximations (RGs, F'*); RGi, r(*'). 



imation, we can still solve this set analytically. We find 



Ef,s,x - - 



UqX 



4(24 + uoA) 



1 



Uo 
24' 



A 



(68) 



The first two terms essentially represent the Hartree term 
whereas the last term originates from the third term on 
the right-hand side of Eq. 111]). As it should be, this 
expression reduces to Egs,\=o — 2 ^^ ^^^ non-interacting 
starting point of the fiow (A = 0). For small Uq, we can 
expand the result for the ground-state energy and obtain 



E, 



1 



1 



1 



gs,A=l 



Un 

2 32 



1536 



ul + Oiul), (69) 



which needs to be compared to the exact perturbative 
results \5U\: 



E\ 



gs 



1 

32 



7 



Mo 



1536 



+ 0{ul). (70) 



Thus, we reproduce the perturbative result at leading or- 
der within the present approximation but not the coeffi- 
cient of the second-order correction. Moreover, it follows 

E„s \=i scales as 



from Eq. (68) that S, 



^gs 



E, 



/Uo 



(71) 



(66) for Uo ^ I. Apparently, this does not agree with the 



result from the large- A^ approximation, see Eq. (54). The 



scaling behavior of the latter has also been confirmed for 
the present model by solving the Schrodinger equation 
numerically, see, e. g., Ref. [5T]. Note that the Hartree 
approximation renders Egg independent of uo in the large 
coupling limit. 

In Fig. |6] we show our results for the ground-state 
energy as a function of the dimensionless coupling uo- 
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For comparison, we also show the exact results [ST] and 
the perturbative results at leading and next-to-leading 
order. Our RG results for the ground-state energy from 
the leading-order approximation is in agreement with the 
leading-order perturbative result at small coupling but 
overestimates significantly the exact results for large val- 
ues of Uq- Thus, the present approximation is insufficient 
and we need to take into account n-point functions of 
higher order, see our discussion below. 

For completeness, we also discuss the behavior of the 
ground-state as obtained from the present approxima- 
tion. First, wc note that p~^' = for A = and n 7^ 0, 
see also Appendix IaI For n = 0, we have P\1,q — > /3/2 
for /3 — > 00. This implies that dxp" = for n ^ 0. 
Thus, only the zero mode p^ contributes to the flow 
of Pgs,A: 



.M 



/3-¥oo 



2ujl + X^ 



(72) 



Note that Pgs,A(''') does not depend on the imaginary 
time T. 

Our results for Pgs,\{T) arc shown in Fig. n\ Since 

(3) 

the three-point function T\ contributes to the flow of 
the ground state pgs already at leading order, we do not 
even recover the leading-order perturbative result with 
our present truncation. For example, this becomes ap- 
parent from the comparison with the exact result in the 
small coupling-limit. The same line of arguments holds 
for the propagator. In fact, not only the three-point func- 

(3) 

tion T\ contributes to the flow of the propagator but 
also the four-point function F 



(4) 



2. Higher-order approximations 

Let us now discuss the role of n-point functions of 
higher order. As discussed in Sect. |llj it can be shown 
that the RG flow of the n-point function depends on the 
flow of the n + 1 and n + 2-point function. Since the 
ground-state is associated with the one-point function p^ 
this means that the two- and the three-point function 
govern the RG flow of pgg. In fact, the most general 
form of the flow equation for the ground state pgs of our 
present model with a constant background potential is 
given by 



(0 



d.p'x 



24 

+ - 



2p-fGl-') 



^A ^A 



3;.o 



k — — c 



G 



W 



(Gf)'(l 



>(3)\ 



{k~kj) 



(73) 



Smce p\ ' = 



for A = and n 7^ 0, only the flow of 

the zero mode p^^ is vanishing and therefore only pj^ 
assumes a finite value. From this observation, it follows 
immediately that the ground state pgs does not depend 
on the imaginary time. This remains true even for d+1- 
dimensional field theories, provided that we consider an 
interaction potential of the form U ^ S{t — r'). In any 
case, from the general flow equation for the ground state 
it is apparent that the three-point function contributes 
to the ground-state at leading order in a perturbative 
expansion. 

Let us now turn to the most general flow equation for 
the propagator (inverse two-point function). This flow 
equation depends on the RG flows of the three- and four- 
point function and can be written as follows 



dxGi"'^ 



Uq 

12 



>2 1 



(^i"')-pE/^i'^i-'(^i"')(f 



/3^ 
1 1 



2 /-/^-vX {-l,rn,-m} I I 
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CO 

Y, J2 G'r"te<"" 



"'^"""■■"'^■'-s^E('Si'"T(«l")(n") 
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A /(4^'T(ff) 
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^E E {G[-^y {Gi'^y ar^^ {r[^^) 



(k./. 



I — — 00 k— — 

00 



3^ E (^A 

k— — oo 



(")^'^Gf)'(fW) 



(k, — k,m, — m) 



^(3) 



[ — l^ — 7n, — k) 



(74) 



r 



Since the right-hand side is proportional to the cou- 
pling Uq, we conclude that we need to include the three- 



The flow equation for the ground state follows from the station- 



ary condition [Sr/5p)\gs = 0. 



and four-point function in our study in order to recover 
the correct perturbative result at leading order. The in- 
clusion of these two correlation functions in the flow of 
the propagator then also guarantees that the results for 
the ground-state energy agree with perturbation theory 
up to second order in the small-coupling limit, see our 
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general discussion in Sect. II C 

For simplicity, we do not include the full flow of the 
three- and four-point function in our numerical studies 
in the following. We rather use a 'static' approximation 
for these correlation functions, i. e. we set their flow 
equations to zero, dx^x — dxT]^ — 0, but we do not 
set the functions themselves to zero. Recall that T\ ' 

and T\^ are not identical to zero at the initial point of 
the RG flow (A = 0), see also Appendix [A| Moreover, we 
set the n-point correlation functions with n > 4 to zero 
as well as their flows. 

Using this approximation, we recover the correct 
leading-order behavior of the propagator and the ground- 
state pgs in the small-coupling limit. Moreover, the 
ground-state energy is correct up to second order in this 
limit, as we have checked analytically. 

In order to go beyond the perturbative small-coupling 
limit, we have to solve the set of flow equations ( 61) , ( 73 1 , 
and ( 74 ) numerically. To this end, we have to trun- 



cate the sum over the Fourier coefhcients, i. e. we 
only take into account the coefficients up to a certain 
value iV,^°"''"^'' . Here and in our subsequent study in- 
cluding an RG improvement, we choose N^™^™ = 100 
and (/3)~^ — T juj — 0.1 for the dimensionless tempera- 
ture. We have also applied this choice to the lowest-order 
RG approximation and compared the results to the cor- 
responding analytic solution. For mq ^ 50, we have found 
that our numerical results for £^gs andpgs deviate from 
the analytic ones by about 1% or lessp^ 

In Fig. [6J we present our numerical results for the 
ground-state energy obtained with this approximation 
(labelled as 'RGs, F^^^^'). The ground-state is shown in 
Fig. [7] We observe that the ground-state energy is now in 
reasonable agreement with the exact result. To be more 
specific, we find that the ground-state energy is less than 
5% larger than the exact value at uq = 10. For uq = 20, 
however, we still have an error of about 10%. In any case. 



.(3) 



.(4) 



the inclusion of V\ and F^^ in our static approximation 
already yields a drastic improvement of our results ob- 
tained with the leading-order approximation. With re- 
spect to the ground-state, we now find good agreement 
between our results and the exact results for uq < 2. For 
larger values of the coupling, we observe that our results 
deviate significantly from the exact results. In fact, also 
the asymptotic functional form of pgg turns out to be 
incorrect. 



3. RG improvement 

Let us now discuss how our RG flows can be improved. 
To this end, we consider the approximation which we 



^* For increasing uq, we have to increase Af^°"''"'' in order to re- 
produce the analytic resuhs. 



have just discussed. The results for ii^gs from this ap- 
proximation are already in reasonable agreement with 
the exact result for uq < 20. For larger values of the 
coupling, this approximation suffers from the fact that 
we have only included F^^ and F^^ 'statically', i. e. we 
have dropped the flow of these correlation functions. Now 
we would like to improve our set of flow equations with- 
out explicitly taking into account the RG flow equations 
for F^^ andF^'*^ 

First, we note that the oscillator frequency to effectively 
changes when the interaction is turned on. To be more 
specific, uj effectively increases with increasing uq due 
to quantum corrections. This becomes apparent when 
we look at the UQ-dependence of the ground state Pgg. 
At A = (non-interacting limit), we have 



Pgs,X=0 



1 

2uJ 



(75) 



for /3 — > oo. Now recall that for increasing A and/or uq, 
the value for the ground state /Ogs,A decreases, see also 
Eq. (72). Thus, a change of the ground state can there- 



fore be viewed as an effective change of the oscillator fre- 
quency. In this spirit, we define an effective A-dependent 
oscillator frequency uj'^' : 



.off. 



1 



(Pgs.xY 



(76) 



In terms of the only non- vanishing Fourier coefficient of 
the ground-state, this can be written as follows: 



-off. 



,cff. 



m 



(77) 



With this effective oscillatory frequency at hand, we can 
improve our set of RG equations for the ground-state and 



(3) 
A=0 



the propagator. To this end we exploit the fact that F 

and ^x=o depend on the oscillatory frequency cu. We now 
replace ut with the effective A-dependent (fiowing) oscilla- 
tory frequency uj'<^' . This can be viewed as an inclusion 
of the three- and four-point function associated with a 
harmonic oscillator with frequency cj'^'- Since ojJ^' is a 
non-trivial function of uq, this improvement allows us to 
include corrections of higher order in a simple manner. 
We would like to add that the change of F)^ and F)^ 
under a variation of A is in principle governed by flow 
equations which we have not taken into account in the 
present study. In this respect, our improvement can also 
be considered as a way to estimate the impact of ne- 
glected n-point functions of higher order. Recall that the 
flow of a given n-point function depends on the n -t- 1- 
and n + 2-point function. This implies that we have in 
general to deal with an infinite tower of flow equations. 
Even if we had included the full flow equations for the 
three- and four-point function, we would still require in- 
formation about the five- and six-point function. Since 
this tower of equations needs to be truncated at some or- 
der, it is always desirable to estimate to some extent the 
effect of the neglected n-point functions. Our proposed 
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RG improvement represents one possibility for such an 
estimate. 

Let us now turn to the (numerical) results from our im- 
proved RG flows. From the construction of these flows, 
it is clear that we still recover the correct perturbative 
result for the ground-state energy up to second order. To 
see this, we note that w^ can be expanded in a power se- 
ries of uq. Our improved three- and four-point functions 
therefore also have a well-defined perturbative expansion: 



A=0 



A=0 



diof- ' 
duo , 



dV 



A=0 



du 



uo + 0{ul) . (78) 



From this expression, it is apparent that the improve- 
ment only affects the third order of the perturbative ex- 
pansion of the RG flow equation for the ground-state en- 
ergy. For the ground state pgs, the second order already 
receives contributions from our RG improvement. The 
latter statement also holds for the propagator. In Fig. [6] 
we show our results for the ground-state energy obtained 
with this RG improvement (labelled as 'RGi, F^")'). The 
results for the ground-state are given in Fig. [7j We ob- 
serve that the improvement brings our RG results for Egs 
and pgs closer to the corresponding exact values. In fact, 
we find that the ground-state energy now deviates less 
than 2% from the exact value at Uq = 10. For uq = 20, 
we have an error of about 5%. We would like to add that 
our RG results for the ground-state energy approach the 
exact results from above when we include n-point func- 
tions of higher order. Comparing our results for pgg with 
the exact values, it even appears that our results are now 
consistent with the asymptotic functional form of pgs for 
large values of the coupling. 



IV. CONCLUSIONS AND OUTLOOK 

We have discussed an RG approach to DFT and an- 
alyzed its properties on very general grounds, including 
its relation to perturbation theory and the computation 
of excited states. A special feature of this approach is 
that it does not rely on a (global) parameterization of 
the density functional. The RG flow rather starts at 
a well-defined analytically accessible starting point. By 
gradually increasing the microscopic interactions, the un- 
derlying RG fiow equation then allows us to follow the 
ground-state of the theory. 

The presented RG approach eventually aims at a study 
of ground-state properties of non-relativistic self-bound 
fermionic systems, such as (heavy) nuclei, from micro- 
scopic interactions. In the present work, however, we 
have restricted ourselves to studies of two simple toy 
models which nevertheless helped us to already test and 
benchmark our approach in a simple but meaningful 
manner for future studies. From a theoretical point of 



view, these two models correspond to a zero and one- 
dimensional field theory. In the present paper, the latter 
was nothing but the quantum anharmonic oscillator. 

Our toy model studies allowed us to illustrate how 
our RG approach works in general and how ground-state 
properties of physical systems can be computed with this 
novel tool. An important feature of our approach is that 
the so-called 2PPI effective action (corresponding to the 
energy density functional) can be computed systemat- 
ically from the underlying microscopic interactions by 
means of a vertex expansion. With respect to nuclear 
physics, this implies that our DFT-RG approach opens 
up a new direction to compute the energy density func- 
tional from, e. g., chiral EFT interactions. Apart from 
the explicit (quantitative) computation of energy func- 
tional and ground-state properties, the presented RG 
approach can be viewed as a tool to gain deeper insights 
into the general structure of density functionals. For ex- 
ample, we have demonstrated how a connection between 
the energy density functional and perturbation theory 
can be established in a simple and systematic fashion. 

For our toy models, we have found that the RG results 
for, e. g., the ground-state energy are in reasonable agree- 
ment with the exact results over a wide range of values 
for the (microscopic) coupling constant, even if we only 
take into account correlation functions up to the four- 
point function. By comparing our results with the exact 
results, we have also demonstrated that the quality of 
the RG results can be systematically improved by taking 
into account n-point functions of higher order. 

The present work should be considered as a starting 
point for various studies. For example, we have al- 
ready mentioned the relation of the present DFT-RG 
approach to the IPI RG approach [IT]. It would also 
be interesting to better understand how it relates to 
other RG approaches widely used to study the nuclear 
many-body problem [55U58) . Apart from these more 
field-theoretically motivated formal studies, another nat- 
ural next step is now to study ground-state properties of 
so-called Alexandrou-Negele nuclei [59] . These are self- 
bound systems in one space and one time dimension, 
consisting of N spinless fermions interacting via a spe- 
cific choice for a long-range attractive and short-range 
repulsive potential U. In this case, the initial condition 
at A = corresponds to a simple harmonic oscillator po- 
tential in which the N lowest lying levels are filled. By 
solving the RG flow equation (11), we then gradually 



remove the background potential V and turn on the two- 
body interaction potential U. The RG results for these 
systems can be benchmarked against those from ab-initio 
MC calculations [5^] and studies with the so-called simi- 
larity RG approach I56j . The present work will help us to 
set up these studies of one-dimensional nuclei |3S] . The 
latter will then help us to further develop and establish 
our novel DFT-RG approach and may pave the way for 
the computation of ground-state properties of realistic 
nuclei from microscopic interactions in a novel alterna- 
tive way. 



16 



ACKNOWLEDGMENTS 

The authors would hke to thank J. Berges, T. Duguet, 
R. J. Furnstahl, H. Gies, F. Karbstein, T. Papenbrock, 
J. M. Pawlowski, J. Polonyi, and A. Schwenk for help- 
ful discussions. Moreover, the authors are very grate- 
ful to R. J. Furnstahl, J. Polonyi, and A. Schwenk 
for useful comments on the manuscript. JB acknowl- 
edges support by the Helmholtz International Center 
for FAIR within the LOEWE program of the State of 
Hesse. Moreover, the authors acknowledge support by 
the Deutsche Forschungsgemeinschaft (DFG) through 
contract SFB 634. 



For the propagator (inverse density-density correlator), 
we find 



where 



Go(r,T')-2Ao(T,T')Ao(T',T), 



oo 



^iw„(T-T') 



/3 



n— — OO 



w,i + w^ 



(A6) 



(A7) 



Note that Ao(r,r') = Aq{t',t). Inserting Eq. ( [A7| ) into 
Eq. ( A5 ) we obtain Fq [pgs] for the non-interacting sys- 
tem: 
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I3uj 
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Appendix A: Initial Conditions 

In this appendix we discuss the computation of the 
initial conditions for our one-dimensional toy model. In 
the spirit of our vertex expansion, we expand the effective 
action about the ground state pgg, see Eq. (14 1. Along 
these lines we also expand the source J about pgs'- 

J{p]- I Ji(r,r')[p(r')~Pgs(r')] 

J t' 

j2(r,r',r")[p(r')-Pgs(T')]x 
x[p(r")-Pgs(r")] + .... (Al) 



Recall that J is a function of r but also a functional of p. 
At the physical ground state pgs, we have J[p = pgs] — 0. 
Using Eqs. (|6| and (14), we then find 



71 - F^') 



J2 -Iq 



(A2) 



Here and in the following we use the index '0' (corre- 
sponding to A = 0) to indicate that we are only dis- 
cussing the non-interacting limit which determines the 
initial conditions for our RG flow equations. 

For a non-interacting theory, the functional Wq is given 
by 



Wo[J]=-\Tv\n(A^^-2J 



1^ 1 



= -oTrlnA- + -E-T^ 



2 ^^ n 



2An- J 



,(A3) 



where Aq denotes the (conventional) propagator and J 
is related to the source J: 



From this expression, we can then extract the ground- 
state energy: 






(A9) 



Next, we compute the Fourier expansion coefficients asso- 



ciated with the propagator G defined in Eq. ( 58 ) . Using 
Eq. ( A6 I , we find 



Gtq - 



— ^^-^cothf^') . (AlO) 

w(a;2-f4a;2) ^ 2 ' 



These coefficients serve as initial conditions for the 
RG flow of the propagator in our study of the one- 
dimensional toy model. 

Let us now turn to the ground state pgs of the non- 
interacting theory. To this end, we consider the relation 
between p and the functional Wq: 

+2 / dT'Ao(r',r)Ao(T,T')J(T') 
Jo 

+ .... (All) 

From this expression, we deduce that 

Pgs(r) = Ao(r, r) = i-coth (^^) . (A12) 

It then follows that the Fourier expansion coefficients de- 
fined in Eq. ( 56 ) are given by 



p£)^Acothf^K„,o. 



2a; 



(A13) 



Moreover, we find 



J{t,t') ^J{t)5{t~t'). 



(A4) 



From a comparison of Eq. ( |A3[ ) with the expansion ( 14 ) 
of the effective action, we immediately obtain Fq: 



robgs] = -TrlnAo^ 



(A5) 



{p{t) - Pgsir)) 

= 2 / rfr'Ao(r', t)Ao(t, t')J{t') + . . . . (A14) 



Plugging the expansion ( Al ) into this expression, we ob- 
tain 
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(p(r) - pg,(r)) - / dr'Go(r,r') 
1 rP rP 



rfr"j-i(T',T")b(r")-Pgs(r")] 



+ 2 j^ dr" J^ dr"'Mr\r'\ r"')b(r") - p,s{r")Mr"') - Pg.(r"')] + 



(A15) 



Here, we have used Eq. (A6) and Go(t, r') — Go{t',t). 
In order to determine the expansion coefRcients j„i we 
compare the left- and right-hand side order by order in 
our functional Taylor expansion in powers of {p — Pgs)- 
For the coefficient ji, we obtain 



dT'Go(T,r')ji(r',r") = <5(T-r"). (A16) 



I 

Thus, we have 

ji(r,r') = Go-i(r,r')=r[,^Hr,r'). (A17) 

Next, we would like to fi.nd a simple relation for F^^ 
This can be done along the lines of our derivation (A17|. 



However, we now need to expand the right-hand side of 
Eq. (A14) up to second order in the source J. From a 



comparison of the expansion coefficients on the left- and 
right-hand side, we then find 



r^^\r,r\r")^Mr,r',r") 

dr'" f dr"" f dr""'Go \r,r"')Ao(r"',T"")Goi(r"",T')Ao(T"",T""')G^i(T""',r")Ao(T""',T"') , (A18) 
Jo Jo 



where we have used Eq. (A17|. The Fourier expansion 



coefficients of the three-point functions can be computed 
straightforwardly and read 
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(A19) 



In our discussion of the one-dimensional toy model, 
we also used the initial condition of the four-point func- 
tion. The computation of this correlation function (as 
well as of correlation functions of higher order) for the 



I 

non-interacting theory can be done following the proce- 
dure detailed above. For the coefficients of the Fourier 
expansion of the four-point function we then obtain 
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where 
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f{uik,ioi.,L^m,Lo) = 640^*" + 48a;'' [So;^ + Atof + AujiLOm + 3w,^„ + 2wfe(2w; + a;„i)] 

+4a;2 [3w^ + 2cjf + 4wfa;„ + lOcjfcj^ + 8cj,a;f„ + 3cj^, + 4w^(2cj; + a;,„) 
+2a;fe(2w( + iOm){Lof + w;a;„i + 2u;,^„) + 2a;^(5u;f + 5w;u;„ + 2a;,^„)] 
+a;fa;^j(w; + a;„)^ + WfeW(W^(u;( + cj„)(2a;j + u;„i) + cj^(wf + w;a;„i + u;,^„) 
+Wfc(2a;; + u;m)(a;f + wjWm + 2[j^) + ^^(wf + 2cjfcj„ + 6a;fu;,^„ + 5a;;w^ + w^) . 



(A21) 
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